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The  major  objective  of  this  project  was  to  develop  automatic  adaptive  grid  methods  for  numeri¬ 
cal  weather  prediction  based  upon  computed  error  estimates.^  We  began  with  the  adaptive  grid 
system  developed  at  Stanford  by  Oliger  and  Berger  for  problems  in  two  space  dimensions  and 
with  a  developmental  hydrostatic  primative  equation  weather  const,  acted  by  Tim  Hogan  at 
NEPRF.  The  adaptive  grid  system  was  extended  to  three  space  dimensions  and  the  data  struc¬ 
tures  of  the  Hogan  model  were  modified  and  a  solver  was  developed  for  use  with  the  adaptive 
grid  system.  Computational  experiments  were  successfully  carried  out  with  this  system  for  a 
barotropic  cyclone  and  an  unstable  baroclinic  jet  with  horizontal  refinement.  The  primitive  hy¬ 
drostatic  equations  were  demonstrated  to  be  inadequate  for  use  with  vertical  refinement  and  the 
equations  derived  by  Browning  and  Kreiss  for  meteorological  applications  were  incorporated 
into  the  model  and  shown  to  reproduce  the  earlier  results  with  horizontal  refinement.  These  re¬ 
sults  are  reponed  on  in  the  paper  by  Oliger,  Skamarock  and  Street  (1989)  and  in  the  Ph.  D. 
Thesis  of  Skamarock  (1986).  We  then  developed  an  interpolation  method  to  preserve  to  near  hy¬ 
drostatic  balance  of  the  atmosphere  for  use  with  vertical  refinement  A  report  on  this  work  is  in 
preparation  and  this  method  is  proposed  to  be  used  in  further  work  for  mesoscale  phenomena 
with  ONR. 


We  have  also  carried  out  several  supporting  investigations  in  the  course  of  this  work.  We  have 
developed  an  accelerated  Schwarz  iteration  for  solving  elliptic  equations  on  adaptive  grid  struc¬ 
tures  for  the  initialization  calculations  needed  in  meteorology.  This  work  is  reported  on  in  the 
paper  of  Oliger,  Skamarock  and  Tang  ( submitted  for  publication )  and  in  the  Ph.  D.  Thesis  of 
Tang  (1987).  We  have  also  investigated  open  boundary  conditions  for  use  with  limited  area 
forecasts.  This  work  is  reported  on  in  the  two  papers  by  Gustafsson  (1987a,  1987b). 
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